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SUMMARY 


Progress in explaining variations of power experienced in the on-line opera- 
tion of a 100-kW experimental wind turbine-generator is reported. Data are 
presented that show the oscillations tend to be characteristic of a wind-driven 
synchronous generator because of low torsional damping in the power train, 
resonances of its large structure, and excitation by unsteady and nonuniform 
wind flow. The report includes a dynamic analysis of the drive-train torsion, 
the generator, passive driveline damping, and active pitch control as well as 
correlation with experimental recordings. Experimental measurements of the 
system transfer function were made by disturbing the blade pitch angle. They 
compared well with the model dynamics up through the frequency of the first 
mode. Oscillations of power were experienced near the first-mode frequency 
and are explained as resulting from the first-mode resonance amplifying distur- 
bances from, for example, blade asymmetries. Control of power about a set 
point used proportional-plus-integral feedback to the pitch actuator, and control 
gains were formulated for reducing disturbances up to frequencies less than the 
first mode. A fluid coupling installed in the high-speed shaft is one solution 
demonstrated for reducing the first-mode resonant amplification. 

A predicted second-mode resonance at 3.5 hertz in the power train was not 
observed experimentally. A small effect of tower motion within the power train 
dynamics (an interaction not modelled) was observed experimentally. Experimen- 
tally observed variations in power at the two-per-rotor- revolution (2P) frequency 
were more than predicted. The larger 2P variation is suspected to be the impact 
of local turbulences not modelled. A wind feed-forward control scheme was em- 
ployed, but failed to attenuate wind- speed-change effects as well as predicted be- 
cause the anemometer measurement of wind speed used in the control was not well 
correlated with the wind speed at the rotor. 


INTRODUCTION 

A 100-kW experimental wind turbine-generator, designated the Mod-0, is 
being used to identify and solve technical problems associated with large wind tur- 
bines. It is located near Sandusky, Ohio and has been operating since November, 
1975 under the federal wind energy program directed by the Department of Energy. 
The National Aeronautics and Space Administration, Lewis Research Center, has 
provided design and test. 

Significant variations of power output are present when the Mod-0 is synchro- 
nized to a utility network (on-line). Some variations occur because of random 
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^ ind speed fluctuations that are too high in frequency for automatic compensation 
through the blade pitch mechanism. Further, blade passage through nonuniform 
winds such as wind shear and tower blockage produces high frequency periodic 
input disturbances. Because the generator assumes a torsion spring-like quality 
when on-line, a spring-mass resonance in the power train can result. This reso- 
nance can be destructive of structure and power quality . 

The oscillation of parallel connected synchronous generators is a persistent 
problem even in current utilities. As discussed in reference 1, problems range 
from power variations of a single generator relative to a utility grid to regional 
oscillations for interconnected networks across the United States. Approaches 
to oscillations from wind turbine synchronous generators are given in references 
2 to 5. The possible implementation of a slip coupling is discussed in reference 2. 
A slip coupling adds damping to the power train but at the cost of more complexity 
and reduced efficiency. In reference 3, transients accompanying wind gusts, with 
a compliant and damped shaft between the turbine and the generator, are calcu- 
lated. A compliant shaft smooths the higher frequency power variations and damp- 
ing prevents resonant oscillations with no steady power loss. However, construc- 
tion of a suitably soft and damped shaft is not easily achieved. The use of a com- 
liant shaft to lower the first resonance to below the turbine rotation frequency is 
analyzed in reference 4; and a method for increased generator damping through 
a variable excitation control is described in reference 5. 

This study examines the nature of the power variations and evaluates, through 
analysis and test, several passive drive-train coupling elements and active pitch 
controls. Only torsional power-train dynamics are considered, and the analysis 
utilizes a linearization of the nonlinear system about an operating point. The wind 
turbine generator is modelled as one machine on an infinite bus with constant gen- 
erator field excitation. The report does not include discussion of induction gener- 
ators, fly wheel inertias, and active damping through variable blade pitch and 
generator excitor controls. 


SYSTEM DESCRIPTION 

An overall view of the Mod-0 wind turbine is shown in figure 1 . It is a two- 
bladed, horizontal axis machine with a rotor diameter of 38.1 meters (125 ft) . 

The axis of the downwind rotor is mechanically yawed normal to the wind. The 
blades are aluminum NACA 23000 series airfoils with a twist of 34° from root to 
tip. The power train components, located inside the nacelle atop the tower, are 
shown in figure 2 . The hub contains the hydraulic servo and bevel gears for col- 
lective pitch control. The low speed shaft supports the weight of the rotor through 
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bearings mounted to the bedplate and is connected to the gear box through a spring- 
type coupling (not shown) . The gear box is a triple mesh of gears that increases 
the 40 rpm rotational speed of the low speed shaft by a factor of 45. The high- 
speed shaft transmits the torque through a parallel set of five V belts, nominally 
having a 1:1 speed ratio. The generator is a four pole, 125 kVA, 480 volt, syn- 
chronous generator rated for 100 kW at a 0.8 power factor. 

Three high speed shaft configurations were analyzed in the test program: a 
stiff steel coupling, a softer spring elastomeric coupling, and a fluid slip coupling. 

A block diagram of the power control feedback loop is shown in figure 3. A 
transducer measurement of the power is compared with set point power and the 
error signal sent to a proportional plus integral structured controller. The con- 
troller commands a pitch angle through the pitch actuator. 

MODEL DESCRIPTION 

The primary goal of the power train modelling was to describe the transfer 
function of generator power-to-blade disturbances for on-line operation. 

Drive Train 

A block diagram of the power train mechanical elements represented as 
lumped masses connected by torsion springs is shown in figure 4. The values of 
stiffness, inertia, internal damping, and external damping used in this model are 
listed in table I. The stiffness and inertia values are from a normal modes analy- 
sis of the Mod-0 wind turbine reported in reference 6. The aerodynamic damping 
of blade inplane motion, parameter Dl, was estimated from curves of rotor power 
at different speeds and is relatively small, especially considering it represents 
the inherent damping available from the wind. The external damping values repre- 
sent viscous losses from friction, windage, etc. which are based on a measured 
80 percent power train efficiency . The relatively larger internal or material 
damping values are estimates based on an assumed 0.05 damping ratio in the 
higher modes. 

Two sets of values are given for the elastomeric shaft. The parameters la- 
beled "stiff" are the parameters for the physical elastomeric shaft tested. The 
set labeled "soft-damped" are for an assumed shaft with an overly optimistic 
damping. Its lower stiffness corresponds to the stiffness obtained from bench 
tests rather than the stiffness measured in actual wind turbine test operation. 

The shaft was instrumented to enable recording its windup while rotating at 1800 
rpm. The higher stiffness recorded is believed to be the result of internal bind- 
ing of the elastomeric shaft steel supports . 
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For the dynamics study, the blade model was a single mass and torsion 
spring. The spring stiffness was sized to give the 2. 17 hertz frequency calcu- 
lated for the blade first inplane mode from structural studies (ref. 7). A nominal 
value of one-half of the rigid body moment of inertia was assigned to the blade tip 
and the remainder to the hub. The sensitivity of this approximation is discussed 
later. 

A block diagram of the generator electrical model connected to an infinite 
bus is shown in figure 5 . The generator model is the voltage-behind-subtransient- 
reactanee model, E", of reference 8. 

Control 

An active generator exciter control is not modelled. Field tests showed the 
exciter response to be slow; hence a constant voltage is used. An active variable 
blade pitch angle control is modelled. The pitch actuator is modelled as a double 
lag at 1.5 hertz based on test measurements of the actuator response. 

The static performance predicted for the blades is given in figure 6. For a 
power-train efficiency of 80 percent, a 125-kW rotor power input produces a 100- 
kW electrical power output. The blade pitch gain, the slope of the power-pitch 
curve, is shown at two wind speeds. The blade pitch gain of 42 kW/deg at the 
higher wind speed is about four times greater than the gain at the lower wind 
speed. To be conservative, that blade gain of 42 kW/deg for the higher wind 
speed is used in controller gain designs. The disturbance gain for wind speed 
changes can also be inferred from figure 6. In the high winds it is 56 kW/m/s 
(26 kW/mph) . In a wind feed forward control scheme to be discussed later, the 
correction in pitch used for a change in wind speed can be obtained from these 
gains as (56 kW/m/s)/(42 kW/deg) = 1.34 degrees of pitch change per meter per 
second change in wind speed. 


Model Equations 

A list of the report symbols, detailed mechanical and electrical model equa- 
tions, and FORTRAN computer program listings are presented in appendixes A, 

B, and C, respectively. 

The power-train model has two nonlinear elements: the stiffness of the low 
speed shaft coupling, and the electrical generator. However, a linear analysis 
was used throughout the study. The low speed shaft was linearized at the oper- 
ating point power level, and the generator equations became linearized by numer- 
ically perturbing the state variables (in both directions) around an operating point. 
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The generator linearization was not particularly sensitive to the perturbation size 
for even a 20 percent change. System transfer functions and frequency responses 
were obtained using the programs of reference 9 . In comparisons of experimental 
and analytical results, a few percent error may exist because the generator power 
and the generator torque were assumed proportional and were used interchange- 
ably. 


MODEL ANALYSIS 

The analysis concerned primarily the transfer function of generator power- 
to-blade disturbances with different plant element dynamics and pitch control. 
Sensitivity studies were made for different blade tip inertias and generator elec- 
trical conditions. 


Plant Dynamics - No Control 

Frequency response magnitude plots of the transfer function of generator 
power-to-blade disturbances are shown in figure 7 . Phase angle plots convey 
mostly redundant information and are not shown. Blade disturbances can result, 
for example, from angle of attack changes due either to changes in wind speed or 
pitch angle. Responses are shown for three high-speed shaft configurations: a 
stiff steel shaft, a soft-damped elastomeric shaft, and a slip coupling set to slip 
2.3 percent at 100 kW. The responses are normalized to one at zero frequency 
and include no pitch control . The stiff steel shaft displays a large peak in the 
response at 0.69 hertz, a frequency very close to the one per rotor revolution 
(IP) frequency of 0.67 hertz. This resonance amplifies wind fluctuations and the 
periodic power variations due to, for example, any differences between blades. 

The peak is less with the slip and elastomeric couplings because of the damping 
the couplings provide. 

From figure 7 it can be seen that the responses of the different configurations 
are about the same at the even multiples of rotor revolution frequency (2P, 4P, 

. . .). The largest difference occurs at the 2P frequency where the slip coupling 
response is about 50 percent lower than that of the stiff shaft. Thus it would be 
expected that the response to disturbances at frequencies above 2P would be nearly 
the same for all configurations. 

The systems' first and second mode natural frequencies and damping ratios 
corresponding to the peaks in the magnitude response plot are listed in table II. 

It is seen that the different couplings affect primarily only the first mode damping 
ratio. The ratio increases from 0.097 for the stiff elastomeric coupling, to 0.10 
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for the stiff steel shaft, to 0.42 for the soft-damped elastomeric coupling, to 0.92 
for the 2.3 percent slip coupling. Note that the damping ratio for the stiff elasto- 
meric shaft is lower than for the stiff steel shaft, showing that a softer coupling 
lowers the resonant frequency but without enough damping may increase the reso- 
nant amplification. 

As discussed, the blade tip moment of inertia was set to one-half of the blade* s 
rigid body moment of inertia. System frequency responses with the slip coupling 
configuration for three blade tip inertias of 0. 1, 0.5, and 1 . 0 of the rigid body in- 
ertia are plotted in figure 8. To keep the results comparable, the blade's natural 
frequency and damping ratio with the hub fixed and the total moment of inertia of 
the blade tip plus hub were held constant. It is seen that the second mode of the 
power train is sensitive to the tip mass. The smaller tip mass results in more 
peaking at a lower frequency. 

In another sensitivity study, changes in the generator and electrical tie line 
conditions were investigated. The sensitivity analysis was for one machine on an 
infinite bus with no active exciter control. The conditions were changed one at a 
time from the power factor of 0. 8 to 1 . 0, the power level of 100 kW to 10 kW, and 
from the external tie line reactance of 0.0103 per unit to 0. 3 per unit. Results 
are considered for the stiff shaft configuration because of its initially low damping. 
The power train's first and second mode characteristics are listed in table 111. 

The first mode damping ratio increases for increased power factor and decreases 
for decreased power. The sensitivity to all changes, however, is small. 

Periodic Response 

The periodic disturbances resulting from linear wind shear and tower blockage 
of the wind were calculated by the rotor aerodynamics program MOSTAB-WT 
(ref. 10). (In the MOSTAB program the blade is divided into 12 segments and al- 
lowed to move in the flap direction. The program assumes a rigid rotating hub 
with no blade torsional dynamics. ) The time domain power variations calculated 
by MOSTAB for three different wind speeds are plotted in figure 9 . In high winds 
the rotor power momentarily drops to less than half. The spectral components of 
the rotor power curves of figure 9 are plotted in figure 10. The curves have con- 
tent only at the even multiples of rotor frequency (2P, 4P, . . . ) which is the case 
for periodic disturbances to two symmetric blades. 

A MOSTAB forcing function and the resulting responses in power output are 
shown in figure 11 . The method used to obtain the power responses in the time 
domain was: 
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(1) Use MOSTAB to calculate the forcing function as a function of time. 

(2) Take a Fourier transform of the forcing function. 

(3) Multiply by the appropriate drive-train transfer function. 

(4) Take an inverse Fourier transform of resulting product. 

The responses obtained for the three high-speed shaft couplings appear similar. 
The output variation displays a 2P and a 4P component nearly equal in amplitude. 
The maximum excursion of power is about 7 kW from the mean. 

The spectral magnitude components of the output power (from step 3 above) 
at the key (2P,4P, . . .) periodic frequencies are plotted in figure 12 for the slip 
coupling configuration with no blade pitch control. These magnitudes are less 
than the MOSTAB disturbance given in figure 10 by the attenuation of the power 
train. Consistent with the time response description, the 4P components are as 
large as the 2P components . 


Closed Loop Pitch Control 

The proportional plus integral blade pitch controller was described in fig- 
ure 3 . System magnitude responses with closed loop pitch control are shown in 
figure 13 for the slip coupling configuration for a family of controller gains each 
having the proportional gain equal to one-half of the integral gain. For a given 
integral gain, that ratio nearly maximizes the damping ratio of the first mode . 
Depending on gain, the control attenuates disturbances to a value between 0.7 and 
0.2 of their open loop value at a frequency of 0.1 hertz. At a frequency a decade 
lower the closed loop amplitude is about a decade lower because of the integral 
control action. All responses peak above the open loop amplitude before reaching 
the system first mode frequency. The peaking develops into a closed loop insta- 
bility with sufficiently high gains at about 0.6 hertz. The effect of the control is 
not significant at the (2P,4P, . . .) blade passing frequencies because the control 
bandwidth is less than the 2P frequency . 

The effect of the controller with the stiff steel shaft is similar to results 
shown for the slip coupling if a filter set to notch out frequencies near the system 
first mode is used in the control. The notch filter used was a second order zero 
with a damping ratio of 0 . 1 (to cancel the system first mode) over a second order 
pole with a damping ratio of 1 . 0 . Of course, the control does not remove the first 
mode resonance, which continues to exist, but the control does not excite this reso- 
nance. Without the notch filter the stiff high-speed shaft coupling is harder to con- 
trol. For example, for system stability there can be no significant proportional 
gain and the integral gain must be low . 
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Wind Feed Forward Control 

A wind feed forward control scheme using the wind speed signal as an open- 
loop input into the pitch angle control was investigaged as a means of attenuating 
the effect of wind speed changes. A schematic of the wind feed forward control is 
shown in figure 14 . The anemometer used was insensitive to the wind direction 
and was located on the nacelle . It was about 4 . 6 meters upwind of the rotor and 
had about a 4. 6-meter of air distance constant. These characteristics are 
modelled with an exponential lead and a lag. Figure 15, a crojs plot of the data 
in figure 14, shows the pitch angle required to keep constant power for different 
wind speeds. The rated and low power curves shown have nearly the same slope, 
1.342 deg/(m/s), thus one gain can approximate all power levels. The feed for- 
ward schedule is cut off for wind speeds below 8 m/ s, and the output of the closed 
loop power control effectively biases the schedule up and down. A potential flaw 
in this scheme is that the wind speed measured by the relatively small anemom- 
eter may not be representative of the rotor average wind velocity. 

The wind speed correction has dynamics because the signal must pass through 
the anemometer and pitch actuator dynamics. With an ideal measurement of the 
wind speed (one that is truly representative of the instantaneous blade average 
wind) and the assumed linear relation between wind speed and power, disturbances 
due to wind speed changes would be attenuated by the following factor 


distance upwind 


' distance constant 
^ wind speed 


wind speed 

\ 2 

s + 1 (actuator constant s + 1) 


The frequency response magnitude of this expression is plotted in figure 16 for 
several wind speeds. It can be seen that the feed forward control theoretically 
should result in greatly reduced low frequency wind disturbance error and error 
reduction up to about 2 hertz. 


RESULTS AND DISCUSSION 

As part of the search for a better understanding of the Mod-0 power variations 
and the means to smooth them, experimental measurements of the power train dy- 
namics were made and compared with model predictions. Measurements were 
made for the three high-speed shaft couplings, a range of controller gains, and 
with wind feed forward control. The transfer function of generator power output- 
to-rotor disturbance was the key relationship sought. 
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Transfer Function Measurement 

To obtain data on the response of output power to rotor disturbances, the 
blade pitch angle was intentionally disturbed. In figure 17, a schematic block 
diagram of the setup used to disturb the blade pitch angle is shown. A signal from 
a random noise generator was summed into the pitch command. The resulting 
pitch disturbance provided data for transfer function measurements. About 
25 minutes of data were recorded to achieve statistically good results for each 
test. The taped signals were either analog or sampled adequately to be nearly 
analog for frequencies of interest below 10 hertz. The reduction of data to its 
spectral content and transfer function relationships was performed using a spec- 
trum analyzer (ref. 11). 

For some tests in days of low winds, the generator was run as a motor driv- 
ing the blades as a fan. Comparisons showed that whether the generator was run 
as a generator or as a motor did not signiiicantly affect the transfer functions 
sought. 

An example measurement of the transfer function of electric power output to 
pitch angle change is shown in figure 18 . The data are most reliable near 0 . 5 and 
2.0 hertz as indicated by the high coherence. That is, the output is attributable to 
the input and not to some other plant noise . Near zero frequency the coherence is 
low and the data are not reliable. The magnitude response does not hold a nearly 
constant value extending to zero frequency because wind speed changes compete 
with the test pitch disturbance. The action of the closed loop control is to sup- 
press power changes. This occurred primarily at low frequencies where the 
wind's content is concentrated (ref. 12), and where the loop 0 ain is highest be- 
cause of the integral control. The coherence is also low near 0.7 and 1.3 hercz 
because of system IP and 2P disturbances not related to the test disturbance. 
Beyond about 2 . 5 cps the coherence is generally near zero, and further results 
are reported to only 2.5 hertz. 

There was no clear measurement of a system second mode resonance which 
the power train dynamic model would predict to be near 3.5 hertz. A second mode 
resonance would be indicated by a peak in the magnitude response accompanied by 
a phase shift of 180° and high coherence. In the sensitivity study, it was noted 
that the second mode peak decreased for increased inertia in the blade tip. How- 
ever, it may not follow that the experimental nonobservance of a second mode can 
be used to increase the tip inertia because the cantilevered and tapered blade was 
represented by only a single mass . 
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Stiff steel shaft , - A transfer function measured in high winds with the stiff 
steel high-speed shaft (no slip coupling) and the corresponding model predictions 
are compared in figure 19. The vertical magnitude axis scale is linear to more 
clearly illustrate the large resonance . The resonant frequency, that frequency 
where the phase is -90° and the magnitude peaks, is seen to be 0.69 hertz which 
is also the calculated value. The experimental resonant peak is lower than that 
predicted by the model. The experimental first mode damping ratio is estimated 
to be 0.22 which is more damping than the calculated value of 0.10 (table II). In- 
creased damping could result fiom slippage in the V belts. About 0.3 percent 
slip was found; but it was not veiy repeatable. A 0. 3 percent viscous slip in the 
model results in a 0.2 damping ratio but lowers the calculated first mode fre- 
quency lower than the measured value. Alternative explanations include presum- 
ing the generator has more damping than calculated. It is more likely, however, 
that nonlinearities become important for large swings in pover. Nonlinearities 
in a system may appear to decrease a resonance because power input at one fre- 
quency may s'lift upon output to different frequencies and thus be missed in the 
data reduction. 

Stiff elastomeric shaft . - A measured and calculated transfer function with 
the stiff elastomeric high-speed shaft is shown in figure 20. A large first mode 
resonant peak, similar to the results with the steel shaft, is evident in the plot. 
The observed resonant xrequency is 0.54 hertz which is also the calculated fre- 
quency. From the transfer function magnitude and phase characteristics, the 
first mode damping ratio is estimated to be 0.22, which is more damped than the 
calculated value of about 0.1. It is noted that experimentally, both the elasto- 
meric and the steel shaft responses give a damping ratio of about 0.22. Slip in 
the V belts in the motoring mode where the elastomeric shaft test was run was 
much l:-ss than the 0.3 percent slip quoted for the steel shaft. The likely explana- 
tion for the greater apparent damping again being that, for large swings in power, 
nonlinearities in the system appear to reduce the resonant peak. 

Slip coupling . - The calculated and measured system transfer functions for 
the slip coupling set to slip 1 . 5 and 4 . 7 percent at 100 kW are compared in fig- 
ures 21 and 22, respectively. The predicted responses are in good agreement 
through the first mode resonance of about 0.5 hertz. Near 2 hertz there exists a 
difference of about 10 dB or a factor of 3. This is discussed next. 

Tower interaction . - In previous transfer function magnitude comparisons the 
amplitude responses of the model were lower than the experimental data near 
2 hertz. Because measurements of tower motion correlate well with power varia- 
tions near that frequency it is believed that there is tower interation . Three ac- 
celerometers were mounted in perpendicular planes on the rear main bearing 
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support of the low speed shaft. Only the horizontal yaw results are reported; the 
other two planes were similar. To show the effect of tower motion on the power 
train, the path is broken into two parts . First, the response of tower acceleration- 
to-pitch is shown in figure 23, and then the response of power-to-tower accelera- 
tion is shown in figure 24. Judging from the coherences, this path is important 
only near 2 hertz where tower natural frequencies are excited. At 1 . 85 hertz it is 
interesting to note that "acceleration-to-pitch" (0.18 g/deg) times "power-to- 
acceleration" (107 kW/g) equals 19 kW/deg. This 19 kW/deg compares well with 
the direct value obtained from figure 21 of about 17 kW/deg measured (and 6 *;W/ 
deg calculated) Thus, with good coherence, the tower motion path accounts for 
the higher response magnitude. Fortunately, tower motion does not cause a seri- 
ous problem for the Mod-0 because it occurs at a relatively high and well placed 
frequency in comparison with the major excitation forces . 

Control gains . - The effect of the controller on the closed loop response of 
power-to-pitch command is shown in figures 25 and 26. Figure 25 was run at 
normal values of proportional and integral gain while in figure 26 the gains were 
increased by a factor of 6. Calculated responses are shown in both cases; it is 
noted that the analytical and experimental results agree well. The response with 
the high gain controller exhibits better low frequency attenuation but it is more 
resonant. The actual system was nearly unstable; it oscillated out of synchronism 
moments after completing the test tape. The model with that high gain was un- 
stable. The phase response calculations are shown to only the frequency of insta- 
bility. Obviously the control gain needs to be some compromise between attenua- 
tion at low frequency and amplification near the first mode. 

Wind feed forward . - Operational data recordings of the slip coupling configu- 
ration were made with and without wind feed forward control. Comparisons varied 
because of wind differences between runs and it was difficult to isolate the effect of 
the feed forward control. A comparison of the coherences of the responses of 
"power-to-the measured wind speed" is shown in figure 27 . The argument for 
wind feed forward is that it. should decrease any relation between wind and power. 

It is seen that in both cases the coherence is mostly low and the transfer function 
data are not reliable. The coherence with feed forward showed a lower relation at 
the very low frequencies but a generally higher relation up to 1 hertz. This would 
indicate that the feed forward as implemented was increasing rather than decreas- 
ing power excursions resulting from unsteady wind conditions over a broad fre- 
quency range. This is attributed to the single, relatively small wind speed sensor 
not being a reliable indication of the instantaneous average wind over the circle of 
area swept out by the blades. 
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Power Time History 

f>ample time history traces of Mod-0 power variations about a 100-kW set 
point are shown in figures 28(a) and (b). The power dropped below the set point 
sometimes because there was not enough wind. In circled areas labeled "over 
power”, the power peaked to near 150 kW. In the figure 28(a), the case without 
wind feed forward, the overshoot was apparently in response to a gust of wind. 

It is noted that control action subsequently reduced the power level although the 
wind speed remained high. In figure 28(b), the case with wind feed forward, the 
power overshoot occurred in response to a drop in the measured wind speed; Ap- 
parently, the feed forward control over-corrected because the measured drop in 
local wind was not representative of the blade average wind speed. Unlike the 
model calculations of figure 11, the power exhibits oscillations at a frequency of 
2P fading in and out The 2P oscillations have at times about a 20-kW zero-to- 
peak amplitude for both low and high winds . This 20-kW (maximum) amplitude 
compares poorly to a predicted (mean) amplitude of only 1.2 kW in a 8.9 m/s 
wind from figure 12 . 


Spectral Analysis of Power 

A frequency domain spectral analysis of generator power measurements is 
shown in figure 29. The peak at the odd IP multiple exists because, as was found, 
the blades were not pitched the same by a difference of 1.7°. The source of the 
peaks labeled "X" at. 2 . 56 hertz and at twice that frequency has not as yet been 
identified. There are peaks at the even multiples of the blade passing frequencies 
(2P, 4P, . . ). The 2P peak differs qualitatively from the others in that it is much 

broader. This may suggest that nonlinear or higher order couplings exist. An- 
other possibility is that the flow field has random local distortions over only part 
of the blade area and lasting only a few revolutions so as to modulate the periodic 
power around the 2P frequency . The buildup around the 2P peak was observed 
"real-time" during the data reduction. It appeared to be random, lending support 
to assumed unsteady local wind turbulence . 

In figure 30, the frequency content in the peaks of the spectrum of figure 29 
are compared with model predictions. The experimental peaks are not a single 
frequency as in the analytical case. The effective content of the broad peaks was 
calculated as the square root of the sum of the squares of individual analyzer 
readings . The analytical points were values from figure 12 divided by the square 
root of two to convert zero-to-peak sinusoidal amplitude to root mean square 
power The comparison shows good agreement at the 4P and 6P frequencies, but 
there is much more content in the experimental data at 2P and 8P. The 8P content 
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was small enough and at a high enough frequency to be insignificant. The average 
2P output power content, about 4 kW rms, is about four times that predicted and 
is discussed further. 

Power oscillation at the 2P frequency has just been shown higher than that 
predicted using input from MOSTAB for a fixed hub applied as a schedule to the 
power train torsional model. An insight is that the 2P power correlates well with 
the blade flatwise bending moment (a nontorsional moment and not in the model) . 
Transfer function data between power and blade flatwise bending moment, from 
experimental data recorded during normal operation are shown in figure 31 . Be- 
cause the coherence is high near the 2P frequency, it cai oe concluded that the 2P 
power variations result probably from a forcing function and not some unknown 
torsional resonance. A 2P forcing function not considered in the MOSTAB calcu- 
lations is nonlinear wind shear. The MOSTAB calculations assumed a linear wind 
shear of plus and minus 15 percent of the hub wind speed. This linear shear 
causes little net aero-torque variations as blade differences cancel. However, a 
quadratic shear with power portional to wind speed would increase the aero-torque 
variations . A calculation that does not include inertial effects from flap dynamics 
shows quadratic shear could about double the 2P spectrum content (add a 7 . 5 kW 
sinusoid to the 2P rotor input or 1 . 3 kW rms to the 2P generator output) . Another 
possible source of 2P power variation not considered in the MOSTAB calculations 
are random local turbulences. 

CONCLUDING REMARKS 

1. Oscillation of electric power near the one-per- rotor- re volution frequency 
(IP) was evident during early on-line tests of the Mod-0 with the stiff high-speed 
shaft. This was found to result from a torsional power-train resonance near the 
IP frequency having only a 0.10 damping ratio. The resonance-amplified IP 
noise resulting from nonsymmetric blade-pitch settings. The control system 
(when without a first-mode notch filter implemented) further increased the oscilla- 
tion . The frequency of the first mode depends upon the power-train torsional 
spring rate . To shift the resonant frequency away from the IP frequency and to 
add damping to the system so as to reduce the amplification, different high-speed 
shaft couplings were studied. The calculated damping ratio of the first mode in- 
creased to 0.42 for a soft-damped elastomeric coupling and to 0.92 for a fluid 
coupling set to slip 2 . 3 percent at 100 kW . Experimental measurements of the 
system transfer function compared well with linear-model dynamics up through 
the frequency of the first mode. In particular, the ability of the slip coupling to 
damp the first mode in comparison with the stiff shaft was demonstrated experi- 
mentally . 



14 


2. Although the model predicted a significant resonance from the second mode 
of the power train, there was no clear measurement proving its existence. 

3. A small tower motion was observed experimentally. However, tower mo- 
tion was not part of the power-train model. 

4. The measured power content near the two-per- revolution (2P) frequency 
was broadband and had some four-times-greater content than predicted. The 
larger 2P power may result from local turbulences acting over only part of the 
blade area. 

5. In a sensitivity analysis, dynamics of the power-train model (with constant 
generator excitation) were not significantly affected by small variations of the gen- 
erator and electrical parameters of power level, power factor, and tie-line reac- 
tance. 

6. Closed-loop pitch controller proportional-plus-integral-feedback gains 
were formulated that could attenuate low-frequency disturbances but only up to 
frequencies less than the system first mode. Excessively high gains drove the 
system toward instability at about 0.6 hertz. 

7 . Analysis of a wind feed-forward control scheme using a wind- speed signal 
to command pitch indicated virtual elimination of low-frequency wind effects and 
attenuation to frequencies well beyond the system first mode. But the large advan- 
tage predicted was not seen, on the average, experimentally because the measured 
wind speed did not represent the instantaneous wind speed past the rotor blades. 
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APPENDIX A 

SYMBOLS 

torque conversion, N-m/P.U. 
exte*" ’al damping, N-m/(rad/sec) 
direct axis 

internal damping, N-m/(rad/sec) 

rmsEMF, P.U. 

field EMF, P.U. 

torque, N-m 

rms current, P.U. 

o 

inertia, N-m/sec 
imaginary number 
constants 

spring rate, N -m/rad 
rotor frequency 
quadrature axis 
resistance, P.U. 
gear ratio 
reference axis 

armature resistance, P.U. 
Laplace variable, sec "* 1 
time constant, sec 
blade torque, N-m 
electrical torque, N-m 
voltage, P.U. 
wind velocity, m/s 
reactance, P.U. 
constant 
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x synchronous reactance, P.U. 

x ac j mutual armature damper circuit reactance, P.U. 

Y admittance, P.U 

p angle between V x and V & , rad 

P pitch angle, deg 

5 angle between quadrature and reference axis, rad 

6 ' angle between E" and reference axis, rad 
6 angle between reference axis and I , rad 

stator rms subtransient flux linkage, P.U. 

0 angle between I a and V_, rad 

a. a 

Subscripts: 
a armature 

d direct axis 

e external 

IC initial condition 

q quadrature axis 

r reference axis 

t terminal 

x conjugate direction to r axis 

0 open circuit 

1 first blade 

2 second blade 

3 hub and low speed shaft 

4 gear box and high-speed shaft 

5 first pulley and "V" belts 

6 second pulley and generator 

«> infinite bus 

Superscripts: 

’ transient 

'' subtransient 

* summed impedance 
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APPENDIX B 
MODEL 

Mechanical Model 

The differential equations of the mechanical model of the Mod-0 power train 
shown in figure 5 are listed below for a perturbation about an operating point. 

*Vl = F 1 " D 1®1 +T b 
J 1^2 ~ F 2 " D 1^2 
J 3®3 = F 3 " F 1 " F 2 " D 3®3 
J /4 = R g F 4 - F 3 - W 
J 5^5 = F 5 " F 4 ‘ D 5®5 
J 6?6 = -F 5 ‘ D 6®6 ‘ [( T e " Te )jo] C 
where the forces between the inertias are 

F l =k l< 0 3-®l> +d l(®3-«l) 

F 2 = k 1^3 “ 6 7 ) + d l(®3 ' ^2) 

F 3 = **3^4 " 0 3* + d 3^4 “ ^3) 

F 4 = k 4 <0 5 ‘ R g®4> + d 4^5 " R g®4) 

F 5 = **5^6 " + d 5^6 " ^5) 


Generator Model 

The differential equations for the generator model shown in the figure 7 block 
diagram are listed next. 
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T 0 = E 'll. + E"I 
e ad q q 


T qO E d “ " E d * < X q - X ") I c 


T do\) “ -X D + E q + (x d " h )J d 
T dO E q = ~ E q (1 + K d> + X xd*d + K d X D + E FD 

S B Vi + K A 

l d = Y[-R(V„ d - E q ) ♦ X"(V„ q - E")] 

I, - -y[x"(V^,- Ejj) * R(V„ q - E")j 

where 

ft = r 0 + R. 

a e 

X’’ - x" + X e 

v «,d " - v « + 6 IC> 

v wq - v„ cos( 20 6 + 6 IC ) 

y - 1 ft ) 2 + (ft”) 2 J 

K 1 = (x " “ ** )/(x d ■ ^ ) 

K 2 " 1 ’ K i 

K d = < x d - x d>< x d - x ^vH " v 2 

X xd = (x d - x J 1 )( x »-x J ,/(x^-x f ) 

Model electrical parameters are given in table IV. The external tie line re- 
sistance and reactance values listed are for the Mod-0 Plum Brook Station location. 
The following four parameters listed were estimated as follows. 
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x" = (x£ + x”)/2 

T dO = W x d> 

T dO = T d( x d /x d ) 

*t =x d‘ x ad 

INITIAL CONDITIONS 

The a. sumed known initial variables were not model states but the machine 
power, the power factor cos <p, and voltage of the infinite bus V . An ini- 
tial condition phase diagram is given in figure 32 . The model values derived 
are in the direct-quadrature axis which is a rotated mirror image of the refer- 
ence bus axis. 

An estimate of current I is 

a 

I s machine power/CV*, cos cp) 

The power into the infinite bus is the machine power minus the transmission line 
ohmic loss. 

2 

bus power = machine power - R I 

6 a 


The current I r in phase with this power is 

I = bus power/V^ 

The derivation of the I component follows. The coordinates of V from 

X a 

figure 32 (or p. 161 of ref. 8) using complex arithmetic are 

V a = V «, + y R e + j X e> 

V a = V » + <*r + iW + j X e> 


V a = (V, + I r R e ' W 




+ L.XJ 
r e' 


The angle /3 between V a and is 
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0 = tan -1 (Imag V a /Real V & ) 


/3 - tan 


-i 'A + 'A 

V oo + h* e - 


The power factor angle <p between I a and V a is the arc cosine of the power 
factor. 

(p - cos - * (power factor) 


From figure 32 


<P - P- 0 


where 0 is the angle of ^ with the reference or 


0 = tan -1 (I x /l r ) 


Thus v 


tan (p = tan(/3 - 0) ^jan ^t a nO 
1 + tan (3 tan 0 


which after substitution for 0 and 13 can be solved for I 


I = -V 

X OO 


% - 4 (- X e * R e «” »l(- 1 r X e + ' r )(V„ * f 


2(-X e + R fi tan <p ) 


The sign of the quadratic used was that for the minimum I x (a negative number) . 
From figure 32 the F" and E qa vectors with respect to the reference axis 

are 

V a * v «, + V R e + i X e> 

E '' = v a + I a< r a + ) x "> 

®qa = ^a + ^ r a + ^ x q^ 
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These variables are related to the model d-q axis. The electrical angle 6 be- 
tween the q and r axes is 


6 = tan _1 (Imag E qa /Real E qa ) 
The angle 6” of the E ,? vector to reference is 

6" = tan~ l (Imag E"/Real E") 


From figure 32 geometry 


Eq - I E"| cos(6 - 5 ") 

Ej = - | E"| sin(6 - 6") 
id = " I I a l sin(6 - 0) 
iq = I i a l cos(6-0) 

The other initial conditions can be found by setting the state equations to zero. 
The Aq state equation and the E q equation are two equations in two unknowns. 

_X D + E q + t x d ~ ^l^d = 0 

with solutions 

E q = E q - K 2< x d - *1 >‘d 
or using the definition of K 2 

Eq * E^ - - x'')I d 


and 


X D =E i + < x d- x |)‘d 
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The state equation solved for Epp is 

E FD = E q (1 + K d) ’ \d l d * K d^D 

The set of numerical values for the figure 32 nominal operating point and pa- 
rameters of table IV are listed below . 


Machine power, P . U 
Power factor .... 

V , P.U 

I a (r,x) 

V a (r,x) 

E qa< r ’ x > 

E”(r,x) 

6, rad 

6'\ rad 

T e * P- U 

I r , P.U 

I X * P U 

p.u 

E^, P.U 

i .ad 


v p - u - 

Eq’ P U - 
V P U. 

E FD , p.u 


«0.8 

0.8 

1 

(0.7953,-0.5898) 
(1.0098,-0.0054) 
. (1.6517,0.8411) 

. (1.1188,0.1225) 

0.4710 

0.1090 

0.8176 

0.7953 

-0.5898 

-0.3985 

1.0525 

0.6381 

-0.8865 

0.4411 

1.0565 

0.9546 

2.8694 
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APPENDIX C 

COMPUTER SUBROUTINES 

Computer subroutines written in FORTRAN generated the initial conditions 
and state derivatives for a linear analysis package. These subroutines are listed 
in figures 33 and 34. 

The FORTRAN symbols which are unlike the report symbols are defined 
below. 


Z(l) - 

Z(11)=0 6 

r * 1 

•CO 

11 

Z(12) = 0g 

03 

CO 

11 

Z(13)=E" 

Z(4) = 0 2 

Z(14) = \ D 

Z(5)=0 3 

Z(15) = E^ 

Z(6)=0 3 

Z(16) = T e (lagged at 10 Hz) 

li 

0 

nT 

Z(17) = p 

Z(8) = 0 4 

Z(18) = p 

Z(9)=0 5 

Z(19) = controller state 

Z{10) =0 5 



CDl=D l , etc. Kll = k^ 

D1 = dp etc. R1 - Rg 

DLT = <5 " TM » T_ 

e IC 

The F vector is the derivative of the Z vector. For example, F(l) = 6^ 
and F(2) = Other constants are 
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B * 2/0.112 sorvo damping 

CONVRT (in-lb/P . U. ) = torque high speed shaft 

Gain (in-lb/P. U. ) = transducer gain ~ V/kW x servo gain (9.4 deg/V) 

150 

x blade gain (42.5 kW/deg) x convert (5869 in-lb/P. U.) 
x speed ratio Rg (45 rpm/rpm 
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TABLE I. - MECHANICAL PARAMETERS 


(a) SI units 


Station 

subscript 

Inertia J, 
2 

N-m-sec 

Stiffness k, 
N-m/rad 

Internal 
damping d, 
N-m 

External 
damping D, 
N-m 

(rad/ sec) 

(rad/sec) 

l 

31 930 

7 682 968 

33 895 

3051 

O 

KJ 

64 686 

a 558 477 , ykW 

2 260 

» 

114 

Four 2.3 per- 

2 859 

0 

131 

.113 

cent slip 





Four steel 

1 154 

75 474 

21 

.113 

Four elastomeric 

1 154 

654 

366 

.113 

(soft-damped) 





Four elastomeric 

1 154 

3 054 

122 

.113 

(stiff) 





Five slip 

1.69 

20 337 

17 

.07 

Five no slip 

.86 

20 337 

17 

.07 

6 

.11 


— 

.18 


(b) U.S 

. customary units 


Station 

Inertia J, 

Stiffness k. 

Internal 

External 

subscript 

lbf-in-sec^ 

lbf-in /rad 

damping d, 

damping D, 




lbf-in 

lbf-in 




(rad/ sec) 

(rad/sec) 

1 

282 600 

68 000 000 

300 000 

27 000 

3 

572 520 

a 4 943 000 VkW 

20 000 

1 010 

Four 2 . 3 per- 

25 300 

0 

1 160 

1.0 

cent slip 





Four steel 

10 210 

668 000 

186 

1 0 

Four elastomeric 

10 210 

5 790 

3 240 

1.0 

(soft -damped) 





Four elastomeric j 

10 210 

27 039 

1 080 

2.0 

(stiff) 





Five slip 

15 

180 000 

150 

.62 

Five no slip 

7.6 

180 OOo 

150 

.62 

6 

24.2 


— 

1 6 


a Linearized about steady state. Cui-ve fit to stiffnesses of low speed shaft 
plus Falk coupling plus gears. 
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TABLE U. - SYSTEM RESONANCES - NO CONTROL 


Configuration 

First mode 

Second mode 


Natural 

frequency, 

cps 

Damping 

ratio 

Natural 

frequency, 

cps 

Damping 

ratio 

2 . 3 Percent slip at 

0.430 

0 920 

3.495 

0.063 

100 kW 





Elastomeric (soft- 

.516 

.420 

3.533 

.058 

damped) 





Elastomeric (stiff) 

.582 

097 

3,522 

Wmm 

Stiff steel 

689 

i 

.104 

3.540 

■a 


TABLE in. - POWER TRAIN (STIFF SHAFT) SENSITIVITY 
TO ELECTRICAL FACTORS 


Variable 

First mode 

Second mode 


Natural 

frequency, 

cps 

Damping 

ratio 

Natural 

frequency, 

cps 

Damping 

ratio 

Nominal a 

0.689 

0 . 104 

3.54 

0.0558 

Power factor = 1 

.713 

.192 

3.55 

.0551 

X e = 03 

.641 


3.50 

.0537 

Power = 11 kW 
(no change in k g ) 

.570 

067 

3.49 

.0684 


a Nominal conditions are power factor - 0 8, external reactance*. 


X D ■- 0.0103, power ■= 103 kW. 

V 
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TABLE IV. - ELECTRICAL PARAMETERS 


sec 


rpff 

V 

, sec 

T d’ 

sec . 

rn M 

L d0 

, sec 

T d0 

, sec 

x d’ 

P.U. 

x d» 

P.U. 

x d’ 

P.U. 

x M 

V 

P.U. 

V 

P.U. 

x ad’ 

, P.U 

x". 

P.U. 

K r 

P.U. 

Kg, 

P.U. 

K d’ 

P.U. 


Xxd* PU 

r a’ P U 

R e* pu 

x e’ P U 

V , P.U 

00 

Base kilovolt-amperes, kVA . . 

Base voltage, V 

Base ohms, ohms per phase wye 


. . 0.0085 

. . . 0.062 

. . . 0.145 

. . . 0.011 

. . . 1.942 

. . . 0.128 

. - . 0.165 

. . . 2.21 

. . . 0.193 

. . . 0.05 

. . . 2.16 

. . 0.1605 

. . . 0.961 

. . . 0.039 

. . . 0.696 

. . . 1.965 

. . . 0.018 

. . 0.00466 

. . 0. 01031 

1 

.... 125 
.... 480 
. . . 1.843 
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Figure Z - Mod-0 turbine generator schematic of nacelle interior. 


Wind, m/s 



Figure 3. * Bfock diagram of closed loop power control. 







Figure 5. - Block diagram of generator third order dynamic model connected to an infinite bus through algebra equations. 











Electrical power 
Rotor power 






Figure 8, - Sensitivity study showing second mode variation in system frequency res- 
ponse magnitude comparisons for three blade tip masses given as fraction of blade 
rigid body moment of inertia. Slip coupling in high speed shaft slip rate of Z 3 per- 
cent at 100 kW. 



Blade angle, deg. 

Figure 9. - MOSTAB rotor power versus blade azimuth 
for one -half revolution. Tower blockage assumed to 
reduce free stream velocity 28 percent over 30-degree 
arc. Blades vertical at 0 P and 181/1 



Figure 1GL - MOSTAB rotor power zero-to-peak sinusoidal spectral components for three wind 
speeds. 



figure li - Time response comparisons of power output with three high speed shaft configurations to 
MOSTAB rotor power input. 



Figure 12. - Calculated spectral content of power output (with Z 3 percent slip at 
100 KW in slip coupling! in response to MOST A B rotor power input for three 
wind speeds. 



Electric power 



Figure 13. - Closed-loop frequency response magnitude comparisons for various pitch 
controller gains K with controller a 027 K [t l/s> r 0. 3}. Plant has slip coupling 
with slip rate of 2. 3 percent at 100 kW. 




Wind, m/s 



| Anemometer | 

I J 

Figure 14 - Block diagram of closed loop power control and wind feed forward control. 
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Figure 15. - Derivation of wind feed forward gain as slope of straight 
line approximation to curves of constant power in pitch angle ver- 
sus wind soeed. 
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Figure 17. - Block diagram of pitch disturbance test for experimental frequency response. 
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Figure 18. - Experimental phase, coherence, and magnitude frequency response 
plot of electric power -to -blade pitch angle, i 05 percent slip at -70 kW. 

0 db * 5 k«V deg. 




Figure 19. - Transfer function of electric power to Wade 
pitch angle for steel high speed shaft +70 kW. 




0 . 



figure 20. - Transfer function of electric power to Wade 
pitch angle for elastomeric high speed shaft -50 KW. 


Experimental 

O Analytical with blade gain 
27. 5 KW/deg 



Figure 2L - Frequency response of electric power to Wade 
pitch angle. 105 percent slip at -70 kW. 0 db - 10 kW/deg. 




Frequency, Hz 

Figure 22. - Frequency response of electric power to blade 
pitch angle. 3,0 percent sli p at -70 kW. 0 db - 10 kW. 






Figure 24. - Transfer function of electric power to horizontal 
acceleration of main bearing. 1. 05 percent slip at -70 kW. 

0 db ■ 750 kW/deg. 




Frequency, Hz 

Figure 25. - Transfer function of electric power to blade 
pitch angle c ommand with nominal control gains of 
proportional - a 025 vfv and integral 1 0. 034 v /vis. 

1. 2 percent slip at 50 kw. 0 dD 1 20 kW/deg. 




Frequency, Hz 


Figure 26. - Transfer function of electric power to Wade 
pitch angle command with high control gains of pro- 
portional * a 137 v/v and integral » a 186 v/v/s. Q 75 
percent slip at +50 KW. 0 db - 20 kVWdeg. 



Wind speed. Power, kW Wind speed. Power, kW 
mphlm/s) mph(m/si 



Frequency, Hz 


Figure 27 - Coherence between generator power and wind speed 
without and with feed forward control Fluid coupling slip of 
2. 3 percent at 100 kW. Proportional control gain of 0. 02 7 v/v 
and integral control gain of 0. 027 sec. 



(a l No feed forward control (10-12-77). 



5CK22. 4M 

lb) With feed forward control UO-8-77). 

Figure 28, - Time history traces of power and measured wind speed with and without 
feed forward control. 2. 3 percent slip at 100-kW. Proportional control gains of 
0, 027 v/v and integral control gam of 0. 034 sec. Over power and 2P oscillations 
noted. 




Frequency, cps 

Figure 29. - 12-12-77 Spectra of generator electric power. 2 3 
percent at 100-kW. Proportional control gains of 0. 027 v/v and 
integral control gain of 0. 034 sec. No wind feed forward con- 
trol. 
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Figure 30. Comparison of analytically predicted 
and measured frequency content of peaks in 
generator power. 






Reference axis, P. U. 


Figure 32. - Phase plot of generator initial conditions for power »0. 8 P. U. , power 
factor • 0. 8. and bus voltage • 1. 0 P. U. I magi nary number j shown as 
used in FORTRAN computer program with complex arithmetic. 


C** MOD-0 INITIAL COMPITIOKS 

SNDRUUTI ME CMSTS(Z, POWER) 

IMPLICIT RLAL(I-L) 

COMPLEX IA,/A,EPP,EUA 
DIMENSION Z(l) 

COItMOil/uEM/Kl, X2, ICO, XXD, XQXPP, XPDXL, Y, R, X, EFD, Til, K3, DELTA, V I MF 
DATA XPPD, XPD, XD, XL/ . 128, . IBS, 2. 210, . 05/, XQ,XPPO/l . 064, . 193/ 
DATA PF, RA/. 8, .013/, RE, XE, /BUS/. 0046O12,. 010307,1./ 

C** GENERATOR CONSTANTS 
100 /CPP=( XPPD+XPPQ) /2 . 

:;pdxl=xpd-xl 

IC1=(XPP-aL)/XP0XL 

K2-1.-K1 

uR=(XD-aPD)*(XPD'aP?)/(XPDXL*XPOXL) 

>;xd=(>:d-xpd)*(x?:‘-xl)/xi dxl 
;:qxpp=xq-xpp 

R=RA*RE 


c** 

200 


c** 


c** 


300 


X=XPP+XE 

Y=1./(P.*R+X*X) 

/ 1 IJF=VBHS 

FI MO IX AS PEP. PC lbl OF A.JDER30M 
TTHEE=TAM(ARCOS( PF) ) 

I R= ( POWEP.-RE* ( ( POWER/ (VIMF*PF) )**2) J/VIMF 


oi|R=SQRT(VlfJF**2-4. *(-XE+RE*TTIlEc)*( 
IXl=(-/IMF*3aR>/(2.*(-XE+Rt*TTHEE)) 

I X2=( - / 1 :IF- oOR)/( 2. *( -XE*RE*TTilEE) ) 
IX=A.ilNl(IXl, 1X2) 

TER* II MAI. CONDITIONS 
I A=CMPLX( lit, I X) 

/ A*V I NF+ 1 A*CI!PLX ( ttt, X F. ) 
iIQA=VA+ I A*Cr!PLX CRA,XQ) 
cPP=VA*l A*C.'iPLX (RA,XPP) 

DELTA=ATAN2 ( A I MAG( tOA) , REAL( i-Oa) ) 

OI.T =ATAN2(AlfiAG(EPP), REAI.(EPP)) 
tPP0=-bA5S(EPP)*SIN(r>lLTA-JLT) 

E PPQ=*CA3S(EPP)*C0S( DELTA- DLT) 
THLTA=ATAM2( IX, IP.) 

I U=-i.ABS( I A) *S I (DELTA- THETA) 
IQ=+CAB5( IA)*COS(DLLTA-T!|ETA) 

GENERATOR. STATES 
EPa=EPPQ-(XPD-:':PP)*|D 
LO=EPQ+(XPO-XL)*IU 
tFO=EPQ*(l. + KO)-;<.Yn*IU-XO*LP 
i f i=EPPQ* I Q+Fpr*p* I D 
Z(13)=E?PD 
Z( 14 ) =LD 
Z( 15 ) =FPO 
Z( IF, )=TM 

LOW SPEED SHAFT STIFFNESS 
1 * 5 = 0 . 04043i-3*S'lRT(A;SS(TM*125. )) 
RETURN 
UW 


ll;*l R*XF+(VINF*PE*I R)*TTHEE*IP) ) 


UKKHNAL PAGE IS 
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Figure 33. - Subroutine CNSTS Fortran listing. 



i.** fKJD-0 JTATE E*l!ATIOMr. 

-llliltOUTiriE ;il »>C(Z^F) 

IflPUCIT REAL(I-K) 

OlilE ISlUfl Z(1),F(J) 

r.on;u):./iiF.'!/n, w., m., rr n, xhxpp, xpuxl, y, r., efu, i m, r.3, ou.ta, •/ 1 mf 

uATA i;n,K!i/.to<,Ei,1.8r5/ # ai,Ca.:7P.T,u,ni Eii«n/‘‘5. # 58G9. # i7.'»,1A./ 
JATA Jl,J3,-4, J5, JG/.2S2GrG, .5/?;».:E6. 25300., 15. ,24.?/ 
bATA CL1, C('3,Lb4, 01'5, CD6/ ,2/E5,1.01E5, 1.,. 52,1.6/ , fiAI?l/7*3391G . / 
JATA .il,U3,D!i,i;:/3.r5,2.E4, lluO. , 150./ . ltZl*)/f»2 .23/ 

LATA TPP t'',7i': , P n ,TPf.0/.0G2, .011,1.042/ 
i. 3TIFF iLiT: J 4*102 14, !(4*6.6oE5,F>4*lSt», J5*7 .57 

51 Fl=i:il*(i(5)-zn))*ul*(z{r 1-2(2)) 

F2=Nll*(Z(5)-Z(i))+L-l*(Z(6)-Z(ft)) 

F3*.C3*(Z( / )-Z( 5) )+l‘3*(±( 0 )-Z( T ) ) 
F4=K4*(z(‘j)-z(7)*:!i)+o4*U(in)-z(n)*ni> 

F5=»;3*(Z(ll)-ZO)) + u5*(Z(12)-Z(10)) 

F(1)=Z(2) 

F(2)*(Fl-C r l*Z(2)) /J1+uAIM*Z( 17)/*ll 
F(3)=Z(4) 

F(4)=(F2-Cri*Z(4))/Jl 

F(5)*Z(6) 

F(G)»(-F1-F2*F3-CP5*Z(G))/J3 

F(7)=Z<3> 

F(3) = (-F3*F4*t‘l-£n4*fa*Z(3))/J4 
F(0)*7.(10) 

F ( 10 ) = (-F4*F5-cr.’5*Z( 10 ) )/J5 
F( 11 )*Z( 12 ) 

ULTA*Z( 11 )*2. +L.ELT-V 

;i*ifp*-\m ir^.i^riTA) 
vi:iFQ*4yif , F*rn5(ULT\) 
tPP0*Kl*Z( 15 )+K2*Z (14) 
lb=Y*(-R*(VI»Fr-Z(15))+X*(VINFQ-Fr > P })) 
ia=-r*(x*(vi,*iFf -z(i3))+ r t*(ynFii-pppi)) 


ORIGINAL PAGE IS 
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rc=EPPO* lvi+Z( 13)* lb 

F(i2)=(-F5*(Ti.-Tr.)*cc;»iyKT-crc*7(i2))/jr> 

F(13)*(-Z(l3)-XnXf>P*lu)/lPP';o 

F( 14)*(-Z( 14)+7( 15)+XPDXL*I i))/l PPL0 

F( 15 )*(-Z( 15)*( 1. +£L))+XXn* I w*tvl)*Z( 14 )+LFu)/TPP0 

F( 1G )=i!Z10*(TE-Z( 1G) ) 

F(17)*Z(18) 

F( 18) «-i)*Z( 1?. )-Ol tEG/w*( Z( 17)-(Z(lii)*0 1 + (Tf !-Z( 1G ) )*BP) ) 
F(10)=T.;-Z(16> 

100 RCTUR'J 
EMD 


Figure 34 - Subroutine HPRC Fortran listing. 
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